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Abstract. The ALMA Observation Support Tool (OST) is an ALMA simulator which is interacted with solely via 
a standard web browser. It is aimed at users who may or may not be experts in interferometry, or those that do not 
wish to familarise themselves with the simulation components of a data reduction package. It has been designed to 
offer full imaging simulation capability for an arbitrary ALMA observation while maintaining the accessibility of 
other online tools such as the ALMA Sensitivity Calculator. Simulation jobs are defined by selecting and entering 
options on a standard web form. The user can specify the standard parameters that would need to be considered 
for an ALMA observation (e.g. pointing direction, frequency set up, duration), and there is also the option to 
upload arbitrary sky models in FITS format. Once submitted, jobs are sequentially processed by a remote server 
running a CASA-based back-end system. The user is notified by email when the job is complete, and directed to 
a standard web page which contains the results of the simulation and a range of downloadable data products. The 
system is currently hosted by the UK ALMA Regional Centre, and can be accessed by directing a web browser to 
http : //almaost . jb.man .ac.uk| 



1. Overview 

The Atacama Large Millimetre/submillimetre Array 
(ALMA0 is an interferometer consisting of 66 dishes cur- 
rently under construction on the Chajnantor plateau of 
northern Chile. Operating between ~90 GHz and ~1 THz, 
it will be the most sensitive instrument in the world at 
these frequencies when completed in 2012. Observations 
with a 16-element ALMA will begin in 2011. 

The Observation Support Tool (OST) provides a new 
method for simulating ALMA images. It is designed to be 
easily accessible to users who may not be interferometry 
experts, and is interacted with solely via a standard web 
browser. No additional software needs to be installed by 
the user, and in contrast to other web-based tools such as 
the ALMA Sensitivity Calculator!! no client-side process- 
ing takes place. Instead, simulation jobs are defined via 
a standard web form and submitted to a remote server. 
The server is running a custom-built script which pro- 
cesses the submitted jobs sequentially, making use of the 
toolkit to perform a full simulation via the gener- 
ation and imaging of a visibility set. When the simulation 
is complete the user receives an email containing a URL 



which points to a web page containing the results of the 
simulation and several downloadable image products. 

The front end makes use of an open source Javascript 
form checking library called LiveValidatiorH which checks 
the simulation parameters in real time. Basic errors can 
thus be trapped before jobs are submitted to the server. 
A second level of server-side error checking is employed 
to test for more complex issues (e.g. sources remaining 
below the horizon at all times) and these error checking 
routines also replicate the Live Validation checks for users 
who choose to disable Javascript in their browsers. 

This article provides a brief overview of the functional- 
ity of the OST by describing the web front-end and the re- 
sults page. The OST is currently hosted by the UK ALMA 
Regional Centre at the University of Manchester Jodrell 
Bank Centre for Astrophysics, and can be accessed by di- 



recting a web browser to http : / / almaost . jb .man. ac .uk. 



1 http://www.almatelescope.org 

http : // almascience . eso . org/ document-and-tools 
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2. The web interface 

The interface is a standard HTML web form with vari- 
ous components of the simulation denned either by text 
boxes or drop-down menus as shown in Figure [T] It is di- 
vided into five main sections, as defined by the leftmost 
column. The column on the right-hand side provides brief 
notes as to the purpose and usage of each item, however 



http : //www. livevalidation. com 
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(evolving) documentation is also available from a hyper- 
link on the web page. The red and green markers show the 
Live Validation library in action, with erroneous parame- 
ters highlighted in red. A brief summary of each section 
on the web interface follows. 

2.1. Array 

The instrument is defined here using a single drop-down 
menu. The available options are two Cycle-0 configura- 
tions (compact and extended), the 12-element Atacama 
Compact Array (ACA) and the full 50-element ALMA. 
Although the full ALMA array will have numerous poten- 
tial configurations there is only one option for it in this 
menu. The reasons for this are explained in Section 12.31 

2.2. Sky Setup 

This section determines the nature of the 'ideal' sky that 
will be used as the input model for the simulation. The 
OST contains a library of example sky models which may 
be used, and flexibility as a simulation tool is provided 
in this section as the user is able to upload an arbitrary 
sky model in FITS format. The only header parameter in 
the FITS file that the user must ensure is accurate is the 
spatial pixel scale. Everything else is ignored by the OST 
and is instead overwritten by the values defined elsewhere 
in the web form, however if the brightness unit keyword is 
Jy/pixcl or Jy/beam then this is also taken into account. 

Right Ascension is conspicuous by its absence next to 
the Declination field, however this was foregone in favour 
of an hour angle parameter (see Section 12. 3[) . 

The final option in this section can be used to re-scale 
the pixel values of the sky model by setting the brightness 
of the peak pixel and scaling all the other pixel values in 
the image array in relation to this. 

2.3. Observation Setup 

The spectral and temporal properties of the observation 
are defined here. The OST at present only does single- 
channel simulations, however the increased sensitivity of- 
fered by large-bandwidth pseudo-continuum observations 
can be simulated by simply entering a large bandwidth, 
together with the central frequency. For spectral line sim- 
ulations the user unfortunately at present has to submit 
each channel as a unique simulation job, however select- 
ing a few representative channels across the model spec- 
tral line cube is suitable for a crude line detection test. 
Frequency cubes can be uploaded, however the OST will 
only make use of the central channel. The user will be 
notified on the results page if this action is taken. 

The required resolution parameter is only taken into 
account by the simulation algorithm if the user selects 
'ALMA' in the Array section. The required resolution is 
used in conjunction with the central frequency to select 
a suitable ALMA configuration from the 28 'out' layouts 



which are bundled with CASA. If the resolution demand 
falls outside the range of what is offered by the most com- 
pact or extended arrays then one of the extreme configu- 
rations is selected and the user is notified. 

A major observing mode for ALMA will be mosaicking. 
In mosaic mode, the OST takes the sky area demands of 
the model into account, and calculates a standard hexag- 
onal mosaic patter The OST will calculate the number 
of pointings required to cover the requested area, but it 
also makes the assumption that the specified on source 
time is per pointing, and not divided amongst them. The 
other option in this section is to simulate a single pointing, 
whereby a pixel-wise attenuation is applied the sky model 
pre-simulation, crudely simulating the primary beam re- 
sponse of the array. The attenuation function is a normal- 
ized Gaussian, the centre of which is placed at the central 
pixel of the sky model. Its full-width at half-maximum is 
an angular value, equal in radians to 

Q PB = 1.224; C 1 ) 
vu 

where c is the speed of light, in metres per second, v is 
the central observing frequency in Hertz and D is the dish 
diameter in metres. 

As mentioned previously, the Right Ascension param- 
eter for the model sky is not present, and instead the user 
must specify a starting hour-angle value. Thus the simu- 
lation is defined in terms of when the source is observed 
in relation to its transit. 

The pointing duration is subsequently specified, and 
there is also the option to specify a number of visits. This 
allows the user to simulate cases whereby a large amount 
of observing time is required but the hour angle ranges 
are stringent. 

Finally the number of polarizations can be selected. 
The OST does not yet perform simulations in full polar- 
ization, and this parameter merely affects the noise in the 
final image. As with spectral lines however, individual po- 
larization maps can be uploaded. This prospect also par- 
tially motivated the allowance of negative brightness val- 
ues in the re-scaling option, although if a negative value 
is entered here then the user is warned on the results page 
in case it was unintentional. 

2.4. Corruption 

Artifacts in an interferometric image derived from a real 
observation originate due to a variety of effects, includ- 
ing calibration errors, atmospheric effects and the thermal 
conditions of the receivers. Even in the case of a perfect 
observation, the latter is something that cannot be mit- 
igated. The thermal noise determines the absolute noise 

5 Pointings in any given row are offset by half the primary 
beam width, and adjacent rows are offset horizontally from 
one-another by one-quarter of the primary beam width, and 
vertically by the primary beam width multipled by sin(60 de- 
grees) 
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Fig. 1. The web interface of the OST, showing the Live Validation library in action. 



'floor' in an interferometric map below which sources can- 
not be detected. 

The RMS of the noise perturbation to the visibil- 
ity, that is the single complex number which is the per- 
polarization, per-channel correlation product of a pair of 
antennas, is given in units of Janskys by: 

a = 2k BTsys (2) 
■qAy/AiJAt 

where ks is the Boltzmann constant, T sys = T rec + T s k y + 
T C mb is the system temperature in K, r\ is the combined 
product of a series of efficiency terms A is the effective area 
of a single antenna in m 2 , Az/ is the channel bandwidth in 
Hz and At is the integration time per visibility in seconds. 

The receiver temperatures have a unique value per 
ALMA band and T cm f, is set to 2.73 K. The sky temper- 
ature is derived from a model of the atmospheric trans- 
parency at the ALMA site via: 



where T atmos is the atmospheric temperature (assumed to 
be 260 K) and 7 is the transmission fraction. 

The single menu option in the 'Corruption' section re- 
lates to three levels of precipitable water vapour (PWV): 
0.5, 1.5 and 2.5 mm. The atmospheric transmission frac- 
tion as a function of frequency for these three levels is 
shown in Figure PJ] The values in this plot have been de- 
rived from the transmission calculator on the Atacama 
Pathfinder Experiment (APEX) web sit^f]. 

Interpolative functions are fitted to these plots so that 
the transparency can be derived for arbitrary frequencies 
and the corresponding sky temperature can be calculated 
via Equation [3] for the selected level of PWV. To account 
for the variation of T s k y within large bandwidths, the val- 
ues are calculated at ten points across the band and an 
average is taken. This value is then added to T rec to form 
T sys for use in Equation [5J 
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Fig. 2. Atmospheric transmission fraction at the ALMA site as a function of frequency for three column densities of 
precipitable water vapour. The vertical bars and the numbers above the plot denote ALMA bands 3-10. 



2.5. Imaging 

The parameters in this section affect how the simulated 
visibilites are Fourier transformed into a sky image which 
then undergoes optional deconvolution. The imaging pro- 
cess for an interferometric data set does not necessarily 
follow a single unique path, thus complete automation of 
this process is a challenge. How the gridded visibilites are 
weighted can affect the final map, with weighting schemes 
offering a trade-off between resolution and sensitivity. The 
OST offers three weighting options: natural weighting, 
pure uniform weighting and an intermediate Briggs (1995) 
weighting scheme. 

When imaging a genuine observation, deconvolution is 
often carried out interactively, allowing the user to adap- 
tively define regions to be targeted by the CLEAN al- 
gorithm and make an educated call as to when deconvo- 
lution should be terminated. If the OST user wishes to 
deconvolve the simulated map then the OST will set a 
termination threshold according to the theoretical noise 
in the map. When the peak of the residual image reaches 
this value then deconvolution ceases. If the deconvolution 
process fails to converge or otherwise fails to reach the 
threshold then it will be terminated by means of a clean 
component limit which is presently set to 10,000. 

The final output image format option in this section 
simply determines whether the OST returns the down- 
loadable data products in FITS format or the CASA im- 
age format. 

3. The results page 

Once the simulation has been processed the OST will send 
the user an email containing a link to a web page similar to 



that shown in Figure [3] The content of this page is influ- 
enced by user feedback and is subject to change at the time 
of writing. The results page is also divided into four sec- 
tions. Sometimes an extra section will appear at the head 
of the page containing warning messages. These generally 
occur when the OST has encountered a non-fatal problem 
with the simulation and has had to take liberties with the 
simulation parameters. One example of this would be the 
requested bandwidth causing the frequency range to spill 
over a band edge, in which case the OST will truncate the 
frequency coverage and notify the user. 

Images are rendered in PNG format and are presented 
with both a linear pixel intensity scale (on the left) and 
with histogram equalization^] (on the right). 



3.1. Overview 

A few useful parameters about the simulation are pre- 
sented here, some of which are repeats of what is entered 
into the webform and some of which are derived (e.g. max- 
imum elevation, resolution of final map, number of point- 
ings). If the simulation used a non-point-source sky model 
then this is rendered and displayed here. 



7 Histogram equalization is an image processing technique 
which adjusts the pixel intensity histogram of an image such 
that it has a flat distribution. This is particularly useful for 
enhancing low level structure in images where the dynamic 
range is governed by the presence of a few very bright pixels. 
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Fig. 3. An example of the results page. 



3.2. CASA simdata 

This section provides three downloads which are designed 
to allow the user to easily transfer their OST simulation 
into the CASA simdata task. The processed sky model is 
offered as a download and the simdata. last file will set up 
simdata with the parameters used for the OST simulation 
by means of the CASA tget command. The pointing file 
contains a list of the directions for each pointing in the 
mosaic. Note however that as mentioned earlier the OST 




Fig. 4. Using the OST on an iPod Touch. Ease of access 
is one of the key features of the simulator. 



uses an hour angle parameter instead of Right Ascension, 
thus the RA of every OST simulation is forced to zero 
hours. 

If the 'single' option is chosen for the pointing strategy 
then in addition to the processed sky model the sky model 
with primary beam attenuation applied is also returned in 
the chosen image format. Such features may be useful to 
users who wish to ignore the simulation components and 
simply use the OST as an online FITS file re-processing 
service. 



3.3. Data products 

The data products section contains solely graphical out- 
put. The most useful image here is probably the render- 
ing of the final simulated map. For mosaics of less than 
30 pointings the pointing directions are overlaid onto this 
image. The link to download the simulated map as either 
a FITS or CASA image is also here. The dirty beam (or 
point-spread function) of the observation is also presented. 
The uv coverage of the simulation is displayed here. This 
is generated by Fourier transforming the PSF rather than 
opening the visibility set and extracting the uvw coor- 
dinates of each measurement. The advantage of this ap- 
proach is that it is much faster, and the colour scale in the 
uv plot gives some indication as to the density of samples 
in a particular region of the uv plane. 
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The frequency set up is distilled into the fourth row of 
plots, showing the frequency range of the simulation in the 
context of ALMA bands 3-10, and the band in which it 
lies. The red atmospheric transmission curve also reflects 
the PWV level that was selected. All three possible trans- 
mission curves are shown in Figure [2] of this document. 

Finally this section also presents a plot of elevation 
against time. Scans are flagged when an antenna drops 
below an elevation angle of 10 degrees. Such scans are dis- 
played in a lighter tone on the elevation plot so the user 
can easily see the fraction of their observation which is 
affected. A message will appear on the results window no- 
tifying the user of any elevation issues, and noise values 
calculated from the on-source time are also scaled accord- 
ingly. 

3.4. Diagnostics 

This section contains only a single result, which is the 
number of seconds that elapse between the simulation job 
being selected from the queue and the job being com- 
pleted. Turn-around times are very favourable, although 
the back-end is being modified to exploit the multiple 
CPU cores of the server, allowing jobs in the queue to 
be processed in parallel rather than the current sequential 
implementation. 
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4. Summary 

The Observation Support Tool is a flexible imaging sim- 
ulator for ALMA which is highly accessible, the software 
requirements it places upon its users consisting of nothing 
more than a standard web browser (e.g. Figure [4]). The 
server-side software which processes the simulations make 
use of the CASA toolkit, and although the OST is not 
a web interface to the CASA simdata task, both systems 
are based on this toolkit and further good agreement is 
ensured via careful matching of assumed parameters. 

The OST itself also serves as an example of the po- 
tential of using remote computing applicatons for radio 
astronomical data processing. The new generation of ra- 
dio interferometers (e.g. EVLA, e-MERLIN, LOFAR and 
eventually MeerKAT, ASKAP and the Square Kilometre 
Array) require high-end computers to perform data re- 
duction, and the use of remotely- accessed HPC facilities 
is likely to become increasingly important, an approach 
which has already been adopted for the processing of 
LOFAR data which takes place on a central cluster fa- 
cility. 

The system is at present in a state of constant refine- 
ment as its performance is checked and user feedback is 
received. The system is completely open for general use 
and has to date processed over 1,500 simulations. 
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